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The topological property of boson-fermion mixture in a one-dimensional optical superlattice is 
studied and the topological insulating phase of interacting boson-fermion mixture characterized by 
a nontrivial Berry phase is identified. The single-particle and boson-fermion exchanging excitation 
spectrums are calculated and we identify the boson-fermion exchanging excitation as the gapless 
topological excitation in the topological phase of the mixture. The different kinds of excitations 
are explained explicitly from the competition among the bulk gap, the on-site boson-boson and 
boson-fermion interactions. The Hamiltonian studied has been partly realized in the state-of-art 
cold atom experiments, and the results are very relevant to the experiments. 
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I. INTRODUCTION 

The discovery of topological insulators has stimulated 
great interests in the studies of topological quantum 
phases m- Besides two- and three- dimensional ones, 
one-dimensional (ID) topological phases also attract in¬ 
tense recent studies [6HT3!. ID systems not only ex¬ 
hibit equally rich topological properties, but also have 
the advantage of reduced complexity. Moreover ID topo¬ 
logical phases are experimentally relevant to the optical 
and photonic superlattices, in which many aspects of the 
topological property have been studied, such as: the lo¬ 
calized boundary states, the adiabatic pumping, the mea¬ 
surement of the Berry phase, et.al. [THED] . 

Besides ID fermion topological phases, ID boson topo¬ 
logical phases have also been predicted in various mod¬ 
els [HHIS]. The studies extend the understanding of 
topological property to systems with different quantum 
statistics. Additionally the underlying models, i.e., Bose- 
Hubbard models, represent one of the simplest systems 
realized in ultracold atom experiments |26[ I27j , based on 
which studying the topological property is an important 
topic currently. Cold atomic systems have also allowed 
the interesting realizations of boson-fermion mixtures, 
which rarely occur in nature. They have been the sub¬ 
ject of considerable experimental and theoretical works 
PSI HH] . It is shown that they exhibit rich physical prop¬ 
erties, including many exotic quantum phases |30H32j . 
Naturally an interesting question rises whether the mix¬ 
tures can possess nontrivial topological properties. 

In the paper, the topological property of boson-fermion 
mixture in a one-dimensional optical superlattice is stud¬ 
ied and the topological insulating phase of boson-fermion 
mixture characterized by a nontrivial Berry phase is iden- 
tihed. For the hardcore-boson case, the topological mix¬ 
ture insulator isn’t affected by the boson-fermion inter¬ 
action. However for the softcore-boson case, the topo¬ 
logical mixture insulator is driven to a trivial one by 
the on-site boson-fermion interaction. We calculate the 
single-particle and boson-fermion exchanging excitation 
spectrums to study the bulk-boundary correspondence. 


In the presence of the on-site boson-fermion interactions, 
the single-particle excitation is gapped. Nevertheless for 
small boson-fermion interactions the single-particle exci¬ 
tation energies are still in the gap and the excited parti¬ 
cle mainly distributes near the boundaries of the system, 
which are the remnants of the topological property. It 
is found that the excitation of boson-fermion exchanging 
is gapless under open boundary condition (OBC) in the 
topological mixture insulator. Moreover the bulk gap in 
such excitation spectrum closes at the topological phase 
transition point. Thus we identify the boson-fermion ex¬ 
changing excitation as the topological excitation in the 
topological phase of the mixture. The two kinds of ex¬ 
citations are explained explicitly from the competition 
among the bulk gap, the on-site boson-boson and boson- 
fermion interactions. Since the Hamiltonian studied has 
been partly realized in the state-of-art cold atom experi¬ 
ments, these results are very relevant to the experiments. 



V 

FIG. 1: (Color online) (a) Schematic illustration of the 

Hamiltonian Eq.Q, which consists of one free fermion chain 
and one interacting bosonic chain coupled by the boson- 
fermion interaction V. U is the strength of the on-site boson- 
boson interaction, (b) The energy gap A and the Berry phase 
7b(f) vs V . 












2 


II. MODEL FOR INTERACTING 
BOSON-FERMION MIXTURE 

We consider the ID Bose-Fermi Hubbard model de¬ 
scribed by the Hamiltonian, 

H = - ^(t,,,+i///,+i+7L.c.) (1) 

i 

— + H.c.) 

i i 

where fi(bi) is the fermionic (bosonic) annihilation and 
creation operators on the i -th lattice site; the hopping 
amplitudes tiy+i are alternating: = ^ 1 (^ 2 ) for odd 

(even) i -th site and are assumed identical for bosons and 
fermions; hf = SJSi and nf = fjfi are the number oper¬ 
ators for bosons and fermions, respectively; U, V are the 
strengths of the on-site boson-boson and boson-fermion 
interactions. Here t 2 = 1 is set as the energy unit. 

In the case of noninteracting pure fermions, the Hamil¬ 
tonian described by Eq.Q is Su-Schrieffer-Heeger (SSH) 
model [33], which writes as, 

HsSH = - + t2flj2i+l + H.C.). 


In the momentum space, it is TLsSHik) = (ti -I- 
t 2 cos k)ax + t2Smkcry. The energy spectrum is E{k) = 
±\/W sinfc)2 -I- (ti +t 2 COsk)^. The gap of the system 
is \ti — t 2 |- The Hamiltonian has the chiral symme¬ 
try <Jz'HssH{k)(Jz = —'Essnik), so its topological in¬ 
variant is a winding number w = — ^9fclnDet(U), 
where V = ti + ^ 26 “*^' is the upper off-diagonal ele¬ 
ment of HssHik). Direct calculations give w = 1 for 
ti < t2 while w = 0 for ti > ^2- Alternately since 
the value of the topological invariant is no more than 
one, its topological property can also be characterized 
by the Berry phase, 7 = ^A{k)dk, with the Berry con¬ 
nection A{k) = i{uk\-^\uk) and Uk the occupied Bloch 
state [331 |3S]. The Berry phase 7 mod 27r takes two 
values: tt for a topological insulator and 0 for a triv¬ 
ial insulator. For the topological phase, there appear a 
pair of zero modes under OBC, which distribute near the 
boundaries (also referred as the boundary states). The 
appearance of the boundary states corresponds to the 
nontrivial bulk topological invariant, which is known as 
the bulk-boundary correspondence. The case of ti = 0 is 
special, when the system breaks into independent dimers 
(two nearest-neighbor sites connected by a bond with the 
amplitude t 2 , see Fig ^ and[^ and it is in the topologi¬ 
cal phase with largest gap. Under OBC, the dimer across 
the boundary is broken into two isolated sites, each of 
which hosts a zero mode. 


In the case of pure bosons, it has been well stud¬ 
ied and the system is in a topological Bose-Mott insu¬ 
lator phase for sufficiently large U at half filling. In 
the hardcore limit U = 00 , the system is exactly solv¬ 
able by mapping the hardcore-bosons to fermions via 
the Jordan-Wigner transformation |3B]. The result¬ 
ing noninteracting fermionic Hamiltonian is the same 

as Eq.|^ of pure fermions, except the term 9\9l = 
—f if L rin^i where the hardcore-bosonic op¬ 
erator 9l = with f'l the mapped 

fermionic operator. Although an additional sign in the 
hopping across the boundary is induced when the num¬ 
ber of hardcore-bosons is even, the topological property 
is not affected. 

In the following, Eq. Q is studied using the exact diag- 
onalization (ED) method under a complete basis of the 
form |{nf,...,nf}„) (g) |{nf,...,nf where nf(nf) is 

the boson (fermion) occupation number of the i -th site 
and a{f3) labels all possible compositions with fixed total 
boson (fermion) occupation number Nb(f) = 
m- Here we consider equal filling factors Nb/L = 
Np/L = 1/2 for fermionic and bosonic species. 


III. THE HARDCORE-BOSON CASE 

The hardcore limit is firstly studied, when no mul¬ 
tiple occupation of bosons is allowed. We focus on 
a nontrivial optical superlattice with ti < t 2 - At 
U = 0 , the fermion and boson subsystems are decou¬ 
pled and the system is a topological mixture insulator. 
Since the hardcore-bosons can be mapped to fermions, 
the Hamiltonian Eq. 0 can be viewed as a fermion 
model with pseudospin. In the large-U limit, the ef¬ 
fective Hamiltonian is the antiferromagnetic Heisenberg 
model: "He// = “ 7 nini+i),where 

ni = rii^i ni ^2 is the total particle number on the i- 

th site, J^i+i = is the nearest-neighbor coupling 

and S, = = 1 , 2 ) is the pseu¬ 

dospin with fiB = fi the fermion operator, fi ^2 = // the 
mapped fermion operator, and a = (a^, CTy, CTz) the Pauli 
matrices. 

To characterize the topological property of the inter¬ 
acting system (V ^ 0), the Berry phase can be defined us¬ 
ing the twisted boundary condition, 7 = i ^{'ijjg\-^\ij;g)d6 
with 0 the twisted boundary phase, which takes value 
from 0 to 27r, and ipg the many-body wave function of 
H (6) [351 [35| . H {9) is the same as the Hamiltonian de¬ 
scribed by Eq.Q except the hopping terms across the 
boundary: /f/i ( 6 f Si), which becomes/f/ie*® (SfSie*®). 
Since there are two chains in a topological mixture in¬ 
sulator, we use -fB(F) = i§{9eB(Fj\d0^\9sB(,F9<^^BiF) 
to characterize the topological property of the mixture, 
where 9b[f) is the twisted boundary phase acquired only 
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FIG. 2: (Color online) (a) The chemical potentials of single¬ 
particle excitations, (b) The excitation energy of hardcore 
boson-fermion exchanging. The inset of (a) shows the distri¬ 
bution of the excited boson. The number of lattice sites is 
L = 8. 
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FIG. 3: Schematic illustration of the conhgurations of the 
ground states with different parameters and fillings under 
OBC in the limit case t\ = 0, where the solid circles (squares) 
represent fermions (bosons). In each conhguration, the up¬ 
per (lower) chain is the fermionic (bosonic) part. The crosses 
represent the sites. Two nearest-neighbor sites connected by 
a bond form a dimer. 


by the bosons (fermions). 

The Berry phase 7 b(f) ^tnd the energy gap of the bulk 
system as a function of the interaction V are calculated. 
The energy gap \s IS. = Ei — Eq with Ei , Eq the eigenen- 
ergies of the first-excited state and the ground state of 
the system under periodic boundary condition (PBC). 
As shown in Figj^(b), the Berry phase remains tt all the 
way, while the energy gap decreases continuously and 
vanishes at 1^ = oo. The results imply that the topologi¬ 
cal invariant of the mixture isn’t affected by the hardcore 
boson-fermion interactions. The nontrivial Berry phase 
characterizing the topological phase usually corresponds 
to gapless topological excitations. Next we calculate 
the single-particle excitation and boson-fermion exchang¬ 
ing excitation to study the bulk-boundary correspon¬ 
dence. The single-particle excitation gap is defined by 
Ac = /r(fV-l-I)—/i(A), where is the chemical poten¬ 

tial computed via ^{N) = Eq{Nb, Np) — EQ{NB — i, Np) 
corresponding to single-boson excitation [or ^i{N) = 
Eo{Nb, Np) — Eo{Nb,Np — 1) corresponding to single¬ 
fermion excitation] with Eq{Nb, Np) the ground state 
energy of iV = Np + Np particles |4T]. The boson- 
fermion exchanging excitation gap is defined by As = 
1/2 — vi with 1/2 = En{L/2 — I, L/2 -|- 1) — Eq{L/2, T/2), 
1^1 = Eo{L/2, L/2) - Eo{L/2 + l,L/2 - 1) [I^. 

We firstly calculate the single-particle excitation en¬ 
ergy. The chemical potentials of adding a particle to 
half—1 and half-filling systems under PBC and OBC are 
shown in Figj^(a) (the excited particle may be hardcore- 
boson or fermion and the results are the same). As long 
as V is turned on, the single-particle excitation is gapped. 
However for small V, the chemical potentials are in the 
bulk gap and the excited particle mainly distributes near 
the boundaries, implying these excitations are closely re¬ 
lated to the boundary states. 

The single-particle excitation spectrum in Fig[^a) can 
be understood in the limit case ti = 0 , when the nonin¬ 
teracting system is deep in the topological phase and the 
boundary state is totally on the boundary site. We con¬ 
sider the case of hardcore-boson excitations (it is sim¬ 
ilar for fermion excitations). At half filling and under 
PBC, each dimer is occupied by one hardcore-boson and 
one fermion with the energy Ei = E — Vv^+i 6 
Appendix, in which the ground-state energies of differ¬ 
ent dimers are explicitly calculated). Since the particles 
behave in a similar way as they do in the SSH model, 
the nontrivial Berry phase remains. Then we take one 
hardcore-boson away and the process can happen in any 
dimer with the excitation energy /r(iV) = Ai -|- I. The 
process of adding one hardcore-boson to the half-filling 
system also happens in any dimer and the excitation en¬ 
ergy is ^j.{N -I- I) = V — 1 — El- 

Next we consider the excitation under OBC and start 
from a system with L/2 — I hardcore-bosons. For H < 3, 
the ground state under OBC is obtained by cutting off 
the dimer with only one fermion, generating the configu- 
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ration shown in Figj^ However for H > 3, it is energet¬ 
ically favorable that the dimer with one hardcore-boson 
and one fermion is cut off. The energy difference of the 
above two processes is: 5E = —1 — Ei, which is negative 
for H >3. Then we add one hardcore-boson to the above 
system and the number of hardcore-bosons becomes L/2. 
For F < 3, it just occupies the boundary state and the 
chemical potential is zero. For F > 3 it instead occu¬ 
pies the dimer with only one fermion and the chemical 
potential is the same as the one under PBC. Finally we 
add one more hardcore-boson. For F < 3, it occupies 
the boundary state where one fermion boundary state 
has already been occupied, thus the chemical potential 
is F. For F > 3, the added hardcore-boson energetically 
favors to occupy the dimer with one hardcore-boson and 
one fermion, and the chemical potentials are the same as 
those under PBC. So the results show that the boundary 
states appear even in the interacting topological phase, 
but their occupations are determined by a competition 
between the topological gap and the interaction. 

Although the single-particle excitations become 
gapped in the presence of the interactions, the excita¬ 
tion of hardcore boson-fermion exchanging is gapless, as 
shown in Figj^(b). Though ’’boson-fermion exchanging” 
can not be performed directly, but it can be effectly re¬ 
alized by the following way: firstly take away a boson 
(fermion) from the system, and then add a fermion (bo¬ 
son) into the system. The results can also be understood 
in the limit case ti = 0. For the PBC case, the changing 
can happen in any dimer and the excitation energies are 
V 2 = —Ei,vi = El. For the OBC case and at half filling, 
the boundary states are occupied by one hardcore-boson 
and one fermion, respectively. If the changing happens 
in the bulk, it is gapped. So it must happen via the 
boundary states, which is gapless and thus is energeti¬ 
cally favored. So the hardcore boson-fermion exchanging 
excitation represents the gapless topological excitation in 
the topological phase of the mixture. In fact it is very 
similar to the gapless spinor excitation at the boundary 
of the ID topological Mott insulator [13]. 


IV. THE SOFTCORE-BOSON CASE 
A. The ED results 

Next we consider the softcore-boson case and the phase 
diagram in the (C/, F) plane is shown in Figj^a). The in¬ 
teraction U drives the boson system from a superfluid to 
a topological Bose-Mott insulator [351 IMl SO] • The 
phase transition happens at very small strength of the in¬ 
teraction U. Then at F = 0 and sufficiently large U, the 
system consists of a fermionic topological insulator and 
a bosonic one, which are independent. So the system is 
in a topological insulating phase of boson-fermion mix¬ 
ture. The further inclusion of the interaction F drives 


the topological phase into a trivial one. We perform the 
calculations at fixed U to show the details. 
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FIG. 4: (Color online) (a) The phase diagram in the ([/, F) 
plane for the softcore case, (b) The energy gap A and the 
Berry phase 7 s(f) as a function of the interaction V. In (b), 
the on-site boson-boson interaction is fixed to 17 = 6. 


The energy gap and Berry phase as a function of F 
are calculated and are shown in Fig|^ (b). A topological 
phase transition occurs at a critical interaction Fc, be¬ 
yond which the energy gap vanishes and the Berry phase 
becomes random. To show the topological property, the 
single-particle excitation energy and the excitation en¬ 
ergy of boson-fermion exchanging are calculated. Similar 
to the hardcore-boson case, the single-particle excitation 
is gapped in the presence of F, and the chemical po¬ 
tentials are in the bulk gap for small F [see Figj^ (a)]. 
However the excitation gaps of boson-fermion exchanging 
under OBC are gapless for F <14, and the topological 
phase transition is accompanied by a closing of bulk gap 
in the excitation spectrum [see Figj^ (b)]. So the boson- 
fermion exchanging excitation represents the topological 
excitation of the topological phase. The results can be 
understood from the competition among the bulk gap, 
the on-site boson-boson and boson-fermion interactions 
in the limit case 4=0. 
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FIG. 5: (Color online) The different kinds of excitation en¬ 
ergies as a function of the interaction V: (a) the single-boson 
excitation; (b) the boson-fermion exchanging excitation. The 
on-site boson-boson interaction is fixed to t/ = 6. The number 
of lattice sites is L = 8. 
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B. Analysis from the limit case under PBC 

Firstly a mixture with JVb = Np = L/2 on a periodic 
chain is considered. For V < Vc-, each dimer is occupied 
by one boson and one fermion, and the energy of each 
dimer is Ei [see Fig[^ (b)]. For V > Vc, the phase sep¬ 
aration is energetically favored and the ground state is 
multi-degenerate. The fermions form trivial Mott insu¬ 
lator, while every two bosons occupy a dimer, with the 
energy E2 = ^ — ■ At the critical point, since ev¬ 

ery two dimers with one fermion and one boson has the 
same energy as the two dimers with bosons and fermions 
separated, we have 2Ei = E 2 . The obtained critical in¬ 
teraction Vc = ^^ 2 E 2 ^'^ — 12.9 for 17 = 6) is consistent 

with the one determined by the Berry phase and the en¬ 
ergy gap shown in Figj^(b). In the topological phase for 
V < Vc, each dimer is occupied by one boson and one 
fermion and the particles hop just like the noninteract¬ 
ing case. So its topological invariant is the same as that 
of the noninteracting case and the system is a topolog¬ 
ical phase of boson-fermion mixture. Then we take one 
boson away. For V < Vc the boson can be taken away 
from any dimer, resulting in a configuration shown in 
the upper one of FigH (a) and the chemical potential is 
fi{N) = Ei-\-l. For V > Vc, if the boson is directly taken 
away from a dimer with two bosons, the resulting dimer 
with one boson has a energy —1. If a fermion in the dimer 
with two fermions is further removed to the above dimer 
with one boson, the energy is lowered to i?i — 1. Thus the 
resulting configuration [the lower one of Fig§ (a)] is the 
ground state. Next we add one boson to the system at 
half filling. For small V the boson can be directly added 
to any dimer, generating a dimer with one fermion and 
two bosons. The energy of such a dimer is E3, so the 
chemical potential is fj.{N + 1) = E 3 — Ei. At a criti¬ 
cal interaction Vd, the first two configurations shown in 
Fig§ (c) have the same energies, i.e. E 3 + Ei = £"4 with 
£4 the energy of the dimer with three bosons, through 
which Vci is determined (14i — 6.87 for U = 6). Then for 
Vci < V < Vc, the second configuration shown in Figj^ 
(c) becomes the ground state and the chemical potential 
is im{N -I- 1) = £4 — 2£i. For V > Vc the boson can be 
added to any dimer with two bosons and the excitation 
energy is £4 — £ 2 . 


C. Analysis from the limit case under OBC 


FIG. 6: (Color online) Schematic illustration of the config¬ 
urations of the ground states with different parameters and 
fillings under PBC in the limit case ti = 0. Except the first 
configuration in (b), all other ground state is multi-degenerate 
and we show only one of them. The number of fermions is 
Np = LI 2 . The symbols are the same as those in Fig[^ The 
red rectangles mark the dimers which the excitation affects 
(see the main text). 


Next we study the ground state and the excitation en¬ 
ergy under OBC. The OBC is obtained by cutting off 
one dimer, but the redistribution of the particles may 
happen to let the system in the ground state. We start 
from the case Np = £/2 — 1. For V <U, the situation 
is the same as the hardcore case [see Figj^a) and Figj^ 
in which the upper two figures have the same configura- 
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FIG. 7: (Color online) Schematic illustration of the config¬ 
urations of the ground states with different parameters and 
fillings under OBC in the limit case t\ = 0. The number of 
fermions is Np = I//2. The symbols are the same as those in 
FigH 


tions]. For U < V < Vc, after cutting off the dimer with 
one fermion and one boson, the boson on the boundary 
site moves to one dimer, while the fermion in the dimer 
moves to the boundary site. The energy cost of such a 
process is SE = E 2 — Ei <0, which is energetically fa¬ 
vored. For V > Vc, the system is in phase separation 
except two dimers in which there are two fermions and 
one boson. The energetically favored configuration is the 
one with one fermion in one dimer and the other two par¬ 
ticles in the other dimer, as shown in Figj^ (a). We then 
add one boson, resulting in a system with Nb = Ll 2 . 
For 17 < 3 ^ the added boson just occupies the bound¬ 
ary state, which corresponds to the zero-energy excita¬ 
tion in the single-particle excitation spectrum shown in 
Fig0b). For 3 < 17 <Vc, the added boson occupy the 
dimer with only one fermion and the chemical potential 
is li(iV) = El 1. For V > Vc, the boson is added to 
the dimer with one fermion and one boson, meanwhile 
the fermion in the dimer moves to the dimer with one 
fermion. The chemical potential is fJ-(JV) = E 2 — Ei + 1 . 
The chemical potentials of the above two cases are the 
same as those under PBC, as shown in Figj^a). Finally 
we further add one boson and the number of bosons be¬ 
comes JVb = L/2 + 1. For small V, the added boson 
occupies one of the boundary state with the excitation 
energy V. From a critical interaction V),' [12], the boson 
begins to occupy the dimer with one fermion and one 
boson. The excitation energy is fj,(N -|- 1) = E 3 — Ei, 
which is the same as the bulk one. For U < V < Vc, 
the boson occupies the dimer with two bosons. The ex¬ 
citation energy is fj,{N -|- 1) = E 4 — E 2 , which is lower 


compared to the bulk one. For V > Vc, the added boson 
enters one dimer with two bosons and the excitation en¬ 
ergy is also fi(N + 1 ) = E 4 — E 2 , which is the bulk one. 
The single-particle excitation energy of fermion can be 
analyzed similarly. Compared to the hardcore case, the 
softcore bosons can condense and the occupation of the 
boundary is also determined by additional competitions 
with on-site boson-boson interaction [/. 

The excitation energy of boson-fermion exchanging at 
half filling can be analyzed similarly. For the case with 
PBC, the changing happens in any dimer [see Fig[^b)] 
and the excitation energies (see the definitions in the pre¬ 
vious section) are: 1/4 = Ei — E2,1^2 = —Ei for V < Vc 
and i^i = —El, V 2 = Ei — E 2 for 17 > Vc. Since 2Ei = E 2 
holds at Vc, the excitation energy has a closing here. 
Next we consider the cases under OBC. For V < U the 
changing happens on the boundary, so the excitation en¬ 
ergies are zero. For V > U the changing happens in the 
bulk. The excitation energy in the region U < V < Vc 
are i^i = Ei — E 2 , 1^2 = Ei — E 2 , which are equal and 
are the same with one branch under PBC. For 17 > Vc, 
vi = —El, V 2 = El — E 2 , which are the same with those 
under PBC. It is notable that ioi U < V < Vc although 
the excitation gap of boson-fermion exchanging is gap¬ 
less, the excitation happens in the bulk. 


V. CONCLUSIONS 


The topological property of boson-fermion mixture in 
a one-dimensional optical superlattice with alternating 
hopping amplitudes is studied. The topological insulat¬ 
ing phase of boson-fermion mixture is identified, which 
is characterized by a nontrivial Berry phase. The single¬ 
particle and boson-fermion exchanging excitation spec- 
trums are calculated and the boson-fermion exchanging 
excitation is identified as the gapless topological excita¬ 
tion of the topological phase of the mixture. The two 
kinds of excitations are explained explicitly from the 
competition among the bulk gap, the on-site boson-boson 
and boson-fermion interactions. The Hamiltonian stud¬ 
ied has been partly realized in the state-of-art cold atom 
experiments and the topological property has been stud¬ 
ied. Besides, the recent techniques of controlling optical 
potential at the single-site level make the studies of the 
boundary states available [M ig. So the results pre¬ 
sented here are very possibly studied experimentally in 
the future. 
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Appendix A: The ground-state energy of the dimers 
with different numbers of particles 


(a) (b) (c) 
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FIG. 8: The dimers with different numbers of particles: (a) 
Nf = 1; (b) Nf = 1,Nb = 1; (c) Nf = 2; (d) Nb = 2; (e) 
Nf = 1,Nb = 2; (f) Nb = 3. The symbols are the same as 
those in Fig® 


We explicitly calculate the ground-state energy of the 
dimers with different numbers of particles. The basis 
we use has the general form \i,j',k,l) with i,j{k,l) the 
fermion (boson) numbers on the two sites of the dimer 
(if the dimer is occupied by only fermions or bosons, the 
basis becomes |i,j) with the fermion or boson num¬ 
bers). 

For the dimer shown in Figj^ (a), the basis is: 
(|1,0), |0,1)), under which the Flamiltonian matrix is, 


Ha 


0 -t 2 A 
-t 2 0 )■ 


Its ground-state energy is —12- The case with one boson 
in the dimer is the same. 

For the dimer shown in Figj^ (b), the basis is: 
(|1,0;0,1),|0,1;1,0),|1,0;1,0),|0,1;0,1)), under which 
the Hamiltonian matrix is. 


/ 0 0 -h -h\ 

0 0 -t2 -t2 

— t 2 —t 2 V 0 

\-t2 -t2 0 V J 


Its ground-state energy is Ei = ^ ^ 


y/v^ + 16tl 
. 2 ■ 

For the dimer shown in Fig® (c), the basis is: (|I, I)). 
The Hamiltonian matrix is 0, so the energy of the dimer 
is 0. If there is one more boson in the dimer, we use the 
basis: (|1,1; 1,0), |1,1; 0,1)) and the Hamiltonian matrix 
writes as. 


Hc = 




Its ground-state energy is V — t 2 . 

For the dimer shown in Fig§ (d), the basis is: 
(|2,0), |1,1), |0, 2)), under which the Hamiltonian matrix 


/ U -V2t2 0 \ 

Hd = I —'/2t2 0 —'/2t2 1 . 

V 0 -V2t2 U ) 

Its ground-state energy is Ei = 

For the dimer shown in Figj^ (d), the basis is: 
(|I,0;2,0),|1,0;I,I),|I,0;0,2),|0,I;2,0),|0,I;1,1),|0,1;0,2)), 
under which the Hamiltonian matrix is. 


2V + U 

— \/2t2 

0 

—t2 

0 

^ \ 

— ^/2t2 

V 

— \/2t2 

0 

—h 

0 

0 

— \/2t2 

u 

0 

0 

— ^2 

— t2 

0 

0 

u 

— \/2t2 

0 

0 

— t2 

0 

— \/2t2 

V 

— '/2t2 

0 

0 

—h 

0 

— y/2t2 

2V + U J 


Its eigenvalues generally can not be expressed in simple 
analytical forms, but can be obtained numerically. Its 
ground-state energy is denoted as E^. 

For the dimer shown in Figj^ (f), the basis is: 
(|3,0), |2,1), |I, 2), |0, 3)), under which the Hamiltonian 
matrix is. 


f 

3U 

— \/3t2 

0 



— V3t2 

U 

— 2^2 

0 


0 

— 2 t 2 

U 

— \/3t2 

V 

0 

0 

— '/3t2 

3U / 


Its eigenvalues are 217 ± a/C/^ -I- 2Uti + — ti,2U ± 

^/U‘^ — 2Uti + Atf+ti and the lowest one corresponding 
to the ground state is E 4 = 2U — + 2Uti + — ti. 
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